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■ Abstract. We investigate spherical accretion to a rotating magnetized 

CTn , star in the "propeller" regime using axisymmetric resistive magnetohy- 

drodynamic simulations. The regime is predicted to occur if the magne- 
tospheric radius is larger than the corotation radius and smaller than the 
light cylinder radius. The simulations show that accreting matter is ex- 
\0 . pelled from the equatorial region of the magnetosphere and that it moves 

away from the star in a supersonic, disk-shaped outflow. At larger radial 
distances the outflow slows down and becomes subsonic. The equatorial 
Q-f matter outflow is initially driven by the centrifugal force, but at larger 

q ■ distances the pressure gradient force becomes significant. We find the 

?-) | fraction of the Bondi accretion rate which accretes to the surface of the 

c/2 . star. 
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1. Introduction 



Rotating magnetized neutron stars pass through different stages in their evo- 
lution (Shapiro & Teukolsky 1983, Lipunov 1992). Initially, a rapidly rotating 
(P < Is) magnetized neutron star is expected to be active as a radio-pulsar. 
The star spins down owing to the wind of magnetic field and relativistic parti- 
cles from the region of the light cylinder tl (Goldreich & Julian 1969). However, 
after the neutron star spins-down sufficiently, the light cylinder radius becomes 
larger than magnetospheric radius r m where the ram pressure of external matter 
equals the magnetic pressure in the neutron star's dipole field. The relativistic 
wind is then suppressed by the inflowing matter (Shvartsman 1970). The ex- 
ternal matter may come from the wind from a binary companion or from the 
interstellar medium for an isolated neutron star. The centrifugal force in the 
equatorial region at r m is much larger than gravitational force if r m is much 
larger than the corotation radius r cor . In this case the incoming matter tends to 
be flung away from the neutron star by its rotating magnetic field. This is the 
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Figure 1. Matter flow in the "propeller" regime for a star rotating 
at f2* = 0.5S1.K* after 6.9 rotation periods of the star. The axes are 
measured in units of the star's radius. The background represents that 
the density and the length of the arrows is proportional to the poloidal 
velocity. The thin solid lines are magnetic field lines. 



so called "propeller" stage of evolution (Davidson & Ostriker 1973, Illarionov & 
Sunyaev 1975). 

The "propeller" stage of evolution, though important, is still not well- 
understood theoretically. We discuss results of axisymmetric, two-dimensional, 
resistive MHD simulations of accretion to a a rotating magnetized star in the 
"propeller" regime. We treat the case when matter accretes spherically with 
the Bondi accretion rate (Bondi 1952). Bondi accretion to a non-rotating and a 
slowly rotating star was investigated by Toropin et al. (Toropin at al. 1999) and 
Toropina et al. (Toropina at al. 2003). Investigation of the accretion to a rotat- 
ing star in the "propeller" regime was started by Romanova et al. (Romanova 
et al. 2003). 
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2. Model 

We simulate the plasma flow in the propeller regime using an axisymmetric, 
resistive MHD code. The code incorporates the methods of local iterations 
and flux-corrected-transport. The code is described in our earlier investigation 
of Bondi accretion to a non-rotating magnetized star (Toropin at al. 1999, 
Toropina at al. 2003). The equations for resistive MHD are 

f| + V-(pv)=0, (1) 

<9v 1 
p- + p(v.V)v = -Vp+-JxB + F, (2) 

f = Vx(vxB, + ^V*B, (3) 

^ + V.(p E v)=-pV-v + -. (4) 
at a 

We assume axisymmetry (d/d(f) = 0), but calculate all three components of 

v and B. The equation of state is taken to be that for an ideal gas, p = 

(7 — l)pe, with specific heat ratio 7 = 7/5. The equations incorporate Ohm's 

law J = cr(E + vxB/c), where a is the electrical conductivity. The associated 

magnetic diffusivity, r/ m = c 2 /(Aira), is considered to be a constant within the 

computational region. In equation (2) the gravitational force, F 9 = —GMpH/R 3 , 

is due to the central star. 

We use a cylindrical, inertial coordinate system (r, (f>, z) with the z— axis 
parallel to the star's dipole moment p and rotation axis fi. The equatorial 
plane is treated as symmetry plane. The vector potential A is calculated so 
that V • B = at all times. The star rotates with angular velocity fi* = fi* z. 
The intrinsic magnetic field of the star is taken to be an aligned dipole, B = 
[3R (p ■ R) — R 2 p]/R 5 with p = p z and vector-potential A = p x R/R 3 . 

We measure length in units of the Bondi radius Rb = GM/c 2 ^, with Coo 
the sound speed at infinity, density in units of the density at infinity p^, and 
magnetic field strength in units of Bq which is the field at the pole of the numer- 
ical star (r = 0, z = Z*). Pressure is measured in units of Bq/8tt. The magnetic 
moment is measured in units of po = BqR? b /2. 

Simulations were done in a cylindrical region (0 < z < Z max , < r < 
Rmax)- The size of the region is less then the sonic radius of the Bondi flow 
Rs'- Rmax = Z max = R s /\/2. Thus matter inflows supersonically to the com- 
putational region. The inflow rate is taken to be the Bondi accretion rate, 
M B = 47rA(GM^) 2 / 9 0O /c^ o , where A = 0.625 for 7 = 7/5. A uniform (r,z) grid 
with 513 x 513 cells was used. 

Initially, the density p(r, z) and the velocity v(r, z) are taken to be the 
values given by the Bondi solution (Bondi 1952) . Thus, initially = 0. Also, 
initially the vector potential A was taken to that of a dipole so that B^ = 0. 
The star was initialized to be rotating at the rate 

At the outer boundaries, r = R ma x and z = Z max , the variables (p, v r , v z , e) 
are fixed and equal to the values for the Bondi flow which has = 0. The in- 
flowing matter is unmagnetized with B = 0. At the inner boundary, the vector- 
potential A = 4> As of the intrinsic dipole field of the star is determined on the 
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Figure 2. Enlarged view of Figure 1, the background represents an- 
gular velocity = v ( f ) (r,z)/r. The bold line represents the Alfven 
surface. Dotted line shows sonic surface. 

surface of the model star. More detailed description of numerical model you can 
sec in Romanova et al. 2003. 



3. Results 

Here, we first discuss simulations for the case = Q*/Qk* = 0.5, (3 = 10~ 7 , 
magnetic diffusivity fj m = 10 . For the mentioned values of (5 and rj m , the 
magnetospheric radius r m Q and corotation radius r cor = (GM/f^) 1 / 3 are equal 
for ~ 0.16. For smaller angular velocities, < 0.16, the matter flow around 
the star is close to that in the non-rotating case. For lu* > 0.16 the flow exhibits 
the features expected in the "propeller" regime. 

Figures 1 shows the general nature of the flow in the propeller regime. Two 
distinct regions separated by a shock wave are observed: One is the external 
region where matter inflows with the Bondi rate and the density and velocity 
agree well with the Bondi solution. The second is the internal region, where 
the flow is strongly influenced by the stellar magnetic field and rotation. The 
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Figure 3. The radial dependence of the different radial forces acting 
in the equatorial plane; F mag = (J x B) r /c is the magnetic force (per 
unit volume), F c = pv^/r is the centrifugal force, Vp = (Vp) r is minus 
the pressure gradient force, F grav = — p(V<& gra v)r is the gravitational 
force, and F tot is the total radial force. 

shock wave, which divides these regions, propagates outward as in the non- 
rotating case Toropin at al. 1999, Toropina at al. 2003). For a rotating star 
in the propeller regime the shock wave has the shape of an ellipsoid flattened 
along the rotation axis of the star. The region of the flow well within the shock 
wave is approximately time-independent. The accretion rate to the star becomes 
constant after about 1 — 2 rotation periods of the star. 

A new regime of matter flow forms inside the expanding shock wave. The 
rapidly rotating magnetosphere expels matter outward in the equatorial region. 
This matter flows radially outward forming a low-density rotating torus. The 
outflowing matter is decelerated when it reaches the shock wave. 

There, the flow changes direction and moves towards the rotation axis of 
the star. However, only a small fraction of this matter accretes to the surface of 
the star. Most of the matter is expelled again in the equatorial direction by the 
rotating field of the star. Thus, most of the matter circulates inside this inner 
region driven by the rapidly rotating magnetosphere. 

The solid line in Figure 2 shows the Alfven surface, where the magnetic 
energy-density equals the thermal plus kinetic energy-density, B 2 /87r = e + 
/cv 2 /2. In the equatorial plane, the Alfven radius is r& ~ 4.5i2*. For r < r^, 
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Figure 4. Angular velocity of matter Q = v^/r versus r in the equa- 
torial region for different angular velocities of the star. The dashed line 
represents the Keplerian angular velocity VLk = ^JGM/r' 6 . 



the magnetic pressure dominates, and the magnetic field is approximately that 
of a dipole. At larger distances the field is stretched by the outflowing plasma. 
Matter inside the magnetopause rotates with the angular velocity of the star. 

Figure 3 shows the radial dependencies of the different radial forces in the 
equatorial plane. One can see that for r > 3.5i?*, the centrifugal force becomes 
dominant in accelerating the matter outward. However, at larger distances, 
r > 5.5i?*, the pressure gradient force become larger and determines acceleration 
of matter. Thus, centrifugal and pressure gradient forces accelerate matter in 
the radial direction. Note, that in most of the region (r > 4.3i?*) the magnetic 
force is negative so that it opposes the matter outflow. 

The results we have shown are for a relatively strong propeller. If the 
rotation rate is smaller, the matter outflow become less intense. For < 
0.16, no equatorial outflow is observed and the shock wave becomes more nearly 
spherical as it is in the non-rotating case. On the other hand the shock wave 
becomes more a more flattened ellipsoid for larger Vt*. Matter rotates rigidly 
inside the Alfven radius r& while at large distances Q, = v^/r decreases, see 
Figure 4. 

Figure 5a shows that the accretion rate to the star decreases as the angu- 
lar velocity of the star increases, M /Mb oc °. Figure 5b shows that the 
accretion rate to the star decreases as the star's magnetic moment increases, 
M I 'Mb oc fi^ 2 ' 1 . We have also done a number of simulation runs for different 
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Figure 5. The left-hand panel (a) shows the fraction of Bondi accre- 
tion rate reaching the star as a function of the star's angular velocity 
u>* = f^/f^if*. The right-hand panel (b) shows the dependence on the 
magnetic moment \x for f2* = 0.5^^. 

magnetic diffusivities in the range f\ m = 1CP 6 — 1CP 4 ' 5 , and from this we conclude 
that M/M B oc {rim) - 7 . 



4. Conclusion 

Axisymmetric magnetohydrodynamic simulations of Bondi accretion to a rotat- 
ing magnetized star in the propeller regime of accretion have shown that: 

(1) A new regime of matter flow forms around a rotating star. Matter falls 
down along the axis, but only a small fraction of the incoming matter accretes to 
the surface of the star. Most of the matter is expelled radially in the equatorial 
plane by the rotating magnetosphere of the star. A low-density torus forms in 
the equatorial region which rotates with velocity significantly larger than the 
radial velocity. Large scale vortices form above and below the equatorial plane. 

(2) The accretion rate to the star is much less than the Bondi accretion 
rate and decreases as the star rotation rate increases (oc Q~ L0 ), (b) as the star's 
magnetic moment increases (oc /U -21 ), and as the magnetic diffusivity decreases 

[« (Vm) - 7 ]- 

(3) Because the accretion rate to the star is less than the Bondi rate, a 
shock wave forms in our simulations and propagates outward. It has the shape 
of an ellipsoid flattened along the rotation axis of the star. 
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